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Ionization of the diffuse gas in galaxies 
Hot low-mass evolved stars at work 
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ABSTRACT 

We revisit the question of the ionization of the diffuse medium in late type galaxies, 
by studying NGC 891, the prototype of edge-on spiral galaxies. The most important 
challenge for the models considered so far was the observed increase of [O iii]/H/3, 
[O ii]/H/3, and [N n]/Ha with increasing distance to the galactic plane. We propose 
a scenario based on the expected population of massive OB stars and hot low-mass 
evolved stars (HOLMES) in this galaxy to explain this observational fact. In the 
framework of this scenario we construct a finely meshed grid of photoionization models. 
For each value of the galactic latitude z we look for the models which simultaneously 
fit the observed values of the [O m]/H/3, [O n]/H/3, and [N n]/Ha ratios. For each value 
of z we find a range of solutions which depends on the value of the oxygen abundance. 
The models which fit the observations indicate a systematic decrease of the electron 
density with increasing z. They become dominated by the HOLMES with increasing 
z only when restricting to solar oxygen abundance models, which argues that the 
metallicity above the galactic plane should be close to solar. They also indicate that 
N/O increases with increasing z. 
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1 INTRODUCTION 



In 1963, iHovle fc Ellis! suggested the existence of an ion- 
ized layer about the Galactic plane, to explain the radio fre- 
quency spectrum observed in directions towards the Galac- 
tic pole. Assuming an electron temperature of 10 4 K, they 
estimated the electron density to be ~ 0.1 cm~ 3 and the 
total mass of ionized gas to be ~ 5 x 10 s Mq. The discovery 
of pulsars (|Hewish et alJll96Sf ) and the use of their disper- 
sion meas ure to estimate th e mean electron density along 
the path (|Rohlfs et all 1 19691 ) indicated the existence of ex- 
tended zones of ionized gas in the Galaxy. The detection of 
faint op tical line emissio n outside the classical Galactic H n 
regions (|Revnolddll97lh clearly established the existence of 
a widespread diffuse ionized medium in the Galaxy. How- 
ever, it was only after the discovery of pulsars in globular 
clusters located far from the Galactic plane that the scale 
height of this war m ionized mediu m (WIM) started being 
really appreciated: iRevnoldsl (1 19891 ) showed that the ionized 
gas extends as far as 4 kpc from the Galactic plane and esti- 
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mated a free-electron scale height of about 1.5 kpc. The de- 
tection of extra-planar i onized gas in edge-on spiral galaxies 
using deep Ha imaging (|Dettmarlll990l : lHoopes et alJll999h 
enabled a better description of the vertical distribution of 
this extraplanar diffuse ionized gas (eDIG) which permeates 
the thick disks and haloes of the galaxies. 

The DIG turns out to be a major comp onent of the 
interstellar medium in galaxies (|Revnolddll991f) and investi- 
gating its topology, composition, ionization conditions and 
relation to the other phases of the interstellar medium - cold 
neutral medium (CNM), warm neutral medium (WNM) and 
hot ionized medium (HIM) - is needed for a complete un- 
derstanding of the evolution of the interstellar medium in 
galaxies and its impact in the intergalactic medium. 

Most specialists agree that massive OB stars in galax- 
ies likely repre sent the main sourc e of ionizing photons for 
the eDIG (see lHaffner et al] |2009i , and references therein). 
However, the fact that the eDIG is detected well above the 
thin disks of galaxies implies some porosity of the interstel- 
lar medium, so that part of the ionizing photons emitted 
by the OB stars, whose scale height is ~ 100 pc or less, 
arrive un absorbed at heights of the order of a kiloparsec 
or more |Hoopes fc Walterbosl [20031 ) . Alternatively, an ad- 
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ditional source of ionization could contribute at such lati- 
tudes. The existence of an additional ionizing source is sug- 
gested by the reported increase of such emission line ratios 
as [N ll]/Ha, [S ll]/Ha, and [O iii]/H/3 with galactic height, 
which cannot be reproduced with models of photoionization 
by hot, massive stars, even taking into account the hard- 
ening of the ionizing radiation due to intervening absorp- 
tion (see lHoopes fc Walterbosll2003l and references therein). 
The sources of additional ioniz ation/heating that are most 



commonly invoked are shocks dCorlins fc Randll200l|). tur - 
bulent mixing layers jSlavin et al .11 19931 ; iBinette et alj2009h . 



magnetic reconnection, cosmic rays, or pho toelectric emis- 
sion from small grains (|Revnolds et al.l l 19991). Note that the 
need of additional heating sources has also be discussed in 
another context, that of ion ized gas in the disks of spir al 
galaxies beyond the H idisk (|Bland-Hawthorn et al.|[l997l ). 

Surprisingly, the role of hot low-mass evolved stars 
(from now on referred to as HOLMES), which are plen- 
tiful in the thick disks and lower haloes of galaxies, and 
have been recognized lon g ago to be a significant source of 
UV photons hi galaxies jHillslll97i iRose fc Wentzellll973l ; 
iTerzianl 1 19741; lLvonlll975|) were reconsidered for the eDIG 
only by ISokolowski fc Bland-Hawthorn! (|l99ll ) but seem to 
have been disregarded since then. However, it has been con- 
vincingly shown that the populations of the HOLMES ex- 
pected to be found in galaxies can e xplain the observed 
line intensities in ea rly-type galaxies (|Binette et al.l Il994l ; 
IStasiriska et al]|2008l ). These stars produce a much harder 
radiation field than massive OB stars, and are therefore able 
to heat the eDIG to higher electron temperatures. 

The purpose of this paper is to show that the expected 
populations of HOLMES, when combined to the popula- 
tion of massive OB stars, are perfectly able to reproduce 
the distribution of emission line ratios in the eDIG. We 
chose to focus on one object: the edge-on spiral galaxy NGC 
891, a galaxy located at a distance of 10 Mp c, with overall 
prop erties similar to those of the Milky Way (jvan der Kruitl 
1 1984 ). that has been extensively observed, especially in opti- 
cal emission lines, providing the most demanding diagnostics 
for our scenario. 

The organization of the paper is as follows. Section [2] 
presents the observational data to be explained. Section [3] 
describes the scenario we explore. Section [4] explains how we 
computed the flux and spectral energy distribution of the ra- 
diation from the massive OB stars and from the HOLMES 
and describes our grid of photoionization models using com- 
binations of those as an input. Section [5] presents the con- 
frontation of our models with the observational data. Finally 
Section [6] summarizes our results and presents prospects for 
future work. The appendix discusses several effects not taken 
into account in our principal grid of models, showing that 
they do not alter our main conclusion, i.e. that the expected 
population of HOLMES can explain the emission lines ob- 
served in the eDIG. 



2 OBSERVATIONAL DATA 

Quite a number of spectroscopic observations of the extra- 
planar gas in NGC 891 have been accumulated over the years 
IIDettmar et al. Hl99d;lDettmar fc Schultdll992i; iPildis et al.l 
1 1994 iRandl 19981 ; lOtte et al.ll200ll ; lRand et al.ll2008f ) as com- 



piled bv lFlores^ Faiardo et aT] (|2009l ) for the DIGEDA0 data 
base. The mo st interesting dat a set from our point of view 
is the one by lOtte et aT] <|200lf ). in which the intensities of 
Ha, H/3, [N 11], [O 11], [O III], and [S 110 were measured along 
a slit placed perpendicularly to the plane of the galaxy. The 
availability of the H/3 line allowed the authors to correct 
the observed line ratios for reddening, assuming an intrinsic 
value of Ha/H/3 of 2.9. Care was taken by the authors to ac- 
count for stellar absorption in the H/3 line as well as possible. 
The m ost important advantage of the data set bv lOtte et al.l 
|200ll) with respect to previous data sets is that it contains 
the [O 11] line, which provides an important diagnostic of 
the ionization conditions in the eDIG, when complemented 
with the [O ill] and hydrogen lines. In absence of [O 11] one 
is forced to rely on either [N 11] or [S 11] to estimate the ex- 
citation level of the gas. However the use of [N 11] implies 
making an assumption on N/O, while nothing is a priori 
known on this abundance ratio within the eDIG. Regarding 
[S 11], while the S/O ratio can be expected to be identical 
to the one in the Sun or in classical H 11 regions, the use of 
this low excitation line in conjunction with [O ill] may lead 
to spurious interpretations: this line is emitted in a zone of 
very low excitation and it is known that, even in classical 
H 11 regions, photoionization models in which the gas distri- 
bution is homogeneous fail to simultaneously reproduce the 
intensities of [O ill], [ O 11] , [N n], [S il] with respect to the 
hydr ogen lines (see e.g. lStasihska et"al 1l2006l ; lLevesque et al.l 
l2010h - 

Figure [T] reports the intensity ratios [O in] /H/3, 
[O n]/H/3, and [N ll]/[0 11] as a function of z, the dis- 
tance to the plane of the galaxy. In absence of t ables giving 
the nu mbers in the original paper, the figures of lOtte et al.l 
(|200ll ) were digitized and the line ratios converted to the 
ones we show in Figure [T] which will be directly compared 
to our models. We also digitized and converted the error 
bars, which are shown in our plots as vertical line33- Note 
that these error bars reflect only th e error s of the flux mea- 
surements. In the original figure of lOttel . there are several 
spikes in [O 11] in the region of the plane of the galaxy. The 
authors warn that it is not clear how physical those peaks 
are. We have therefore removed the corresponding points 
in our figure. As has been noted by several authors, the 
[O m]/H/3 ratio increases steadily with \z\. The [O n]/H/3 
ratio also shows this tendency, although less strongly. The 
[N n]/[0 11] shows a shape in "M", strongly suggesting that 



1 DIGEDA (Diffuse Ionized Gas Emission Data base) is available 
in Vizier website: |http:/ / vizier. u- strasbg,fr/cgi-b in/VizieR| 

2 Unless explicitly stated otherwise, in the whole paper [N n], 
[O 1], [O n], [O III], and [S 11] stand for [N 11] A6584, [O 1] A6300, 
[O 11] A3726 + A3729, [O in] A5007, and [S 11] A6716, respectively. 

3 T he error bar s were obtained from the errors plotted in Fig. 
6a of lOtte! j200lT) by assuming that the errors in [O Ill]/H,3 and 
[O n]/H/3 were 2.9 times those in [O nrj/Hc? and [O Il]/Ha re- 
spectively, and that those in [N 11] /[O 11] were the quadratic sum 
of the error in [N 11] /Ho an d oii/Hq. 

4 The digitized data from lOtte et"aL I j200ll) and the data for 
the error bars were added to an updated version (V2.0) of 
DIGEDA. In this new version we corrected a confusion be- 
tween [S n]A6717 + 6731 and [S ll]A6717 and added a code, 
OBS_ID, for the identification of each error with each observa- 
tion (OBS.ID crror =OBSJD data + 100000). 
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Figure 1. Observed values of [O m]/H/3, [O ii]/H/3, and [N ll]/[0 ll]in NGC 891, corrected for reddening, as a function of the distance 
z to the galactic plane. The data are from lOtte et al. I fcooily 



the N/O ratio is not constant. Ignoring this fact when in- 
terpreting the data in terms of the ionization state of the 
DIG would then lead to spurious results. Note that part of 
the light from H uregions in the disk could be dust-scattered 
outside the disk and contribute to the extra planar emission. 
Since the albedo is almost independent of wavelength in the 
optical the line ratios of the scattered component should be 
similar to the ones near the galactic plane. The fact that the 
extra planar emission has different line ratios implies that 
the scattered com ponent is not dominant. 

iRandl (|l998l ) also observed the faint [O i] l ine (with a 
slit n ot exactly at the same position as that of lOtte et al.l 
l200ll ) and plotted it as a function of z. Since the [O 1] line 
has a similar problem as [S 11] when attempting to interpret 
observations with the help of homogeneous photoionization 
models, we will not use this observation to constrain our 
models. However, once our fitting procedure is achieved, we 
will compare our predictions with those observations (see 
Sect [5]). The same will be done with the He I A5876 line, 
extr acted from Fig. 3 of iRand et all (120081). 

iRand et all (|2008l ) and lRand et all (|201lh reported mid 
infra-red spectroscopy of the diffuse ionized halo in NGC 
891, which showed that the [Ne ill] A15.6/xm/[Ne n] A12.8/zm 
is enhanced in the extraplanar pointings relative to the disk 
pointing. We will also confront our models to these observa- 
tions (see Sect [5]). 

Finally, observations of the verti cal di st ributi on of 
the Ha surface bright ness given by IRandl (jl998T ) and 
iMiller fc Veilleuxi (|2003h provide important information on 
the number of ionizing photons absorbed as a function of 
height above the galactic plane. 



3 SCENARIO AND MODELLING POLICY 

We first present a qualitative view of the scenario that we 
propose to explain the emission line features of the DIG. 
This is shown in Fig. [2] where the galactic plane is located at 
the bottom of the figure. The extraplanar gas is distributed 
in clouds, between which ionizing photons emitted by the 
massive stars located in the disk of the galaxy can travel up 
to high latitudes. For simplicity, we represent the gas clouds 
by rectangles. Ionized gas is represented in blue while neu- 
tral gas is represented in red. For both components, a darker 




Figure 2. Schematic representation of the extraplanar gas and 
ionizing stars. Ionized gas is represented in blue, neutral gas in 
hashed-red, with a lighter shade for less dense regions. See text 
for a detailed description. 



shade indicates higher density gas. All that is white repre- 
sents very diffuse, coronal gas which is neither emitting in 
the optical nor absorbing. The big blue stars in the plot 
represent massive OB stars. They are located close to the 
galactic plane, and surrounded by H n regions, some of which 
are density-bounded and allow leakage of ionizing photons. 
The small purple stars represent the HOLMES, which can 
be either central stars of present-day planetary nebulae, or 
hot pre-white dwarfs whose evolution time-scale is larger 
than that of the nebular envelope they ejected in a previous 
stage. These stars are much hotter than massive OB stars. 
They are also much more numerous and their height-scale is 
that of the old stellar population in the galaxy. Most plan- 
etary nebulae, although surrounded by a gaseous envelope 
of density n c ~ 10 2 — 10 4 cm -3 are density-bounded and al- 
low leakage of hard ionizing photons in the diffuse medium; 
some o f them, however, are ionization-bounded (see e.g. Fig 
3 from iMarigo et al. | |2004 which shows the expected pro- 
portion of optically thick and thin planetary nebulae among 
a simulated sample of bright planetary nebulae). In both 
cases, the ionizing photons produced by the central stars 
are ultimately absorbed by gas - be it the proper planetary 
nebula or the diffuse interstellar gas - and, at the spatial 
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resolution of the observations, they can be counted as pho- 
tons available to ionize the eDIG. More importantly, most 
of the HOLMES are not surrounded by a planetary nebula, 
and the total number of ionizing photons they produce is 
thus far greater than what can be estimated from planetary 
nebulae alone. 

The part of the radiation from the massive stars em- 
bedded in the classical H uregions escapes and, in a first ap- 
proximation, ionizes the "bottom" parts of the gas clouds. 
As z increases, the surface fraction of the gas clouds that are 
not shielded by intervening neutral zones becomes smaller. 
At the same time, the radiation originating from the massive 
stars becomes harder, since it is the photons with energies 
closest to the ionization threshold of hydrogen which are 
absorbed first. 

The HOLMES, which are distributed in the galaxy thick 
disk and halo, also contribute to the ionization of the clouds. 
Their influence with respect to that of OB stars increases 
away from the galactic plane, and the "bottom" of the 
clouds is ionized by a radiation field that becomes harder 
and harder. At any latitude, HOLMES will also produce an 
ionized skin on all the sides of the gaseous clouds that are 
not facing the OB stars. The clouds have to be optically 
thick, i.e. to possess a neutral core, since the [O i] line is 
detected at all heights above the galactic plane. 

In the photoionization models presented below we con- 
sider, for simplicity, that each cloud is ionized only from 
the "bottom", by the summed radiation field from the OB 
stars and from the HOLMES|f|. The models are computed in 
the plane parallel approximation, assuming that each cloud 
is of constant density. The computations are stopped when 
the ionization fraction of hydrogen drops below 0.1%. Fig. [3] 
shows the details of the ionization structure of a given cloud. 

Note that the covering factor of the whole population of 
clouds is not necessarily unity, and some ionizing radiation 
may escape from the galaxy. 



4 QUANTITATIVE MODELLING 

In this section, we first explain how we compute the ionizing 
radiation field that reaches the clouds, then we describe the 
grid of photoionization models constructed to interpret the 
observations. 



4.1 OB stars 

The ionizing flux issuing from OB stars is obtained from the 
evolutionary stellar p opulation synthesis code Starburst99 



|Leitherer et al. 19991). W e adop ted the stellar initial mass 
function of Kroupa et alj (1 19931 ) for stellar masses ranging 



from 0.1 to 100 Mq . The stellar tracks are the Geneva ones 
|Maeder fc Mevnet|l994| ) with enha nced mass loss. T h e stel - 
lar model atmospher e s are those of lPauldrach et al] (12001 ) 
and iHillier fc Milled (|l998l ), implemented by ISmith et all 
|2ooa )" in Starburst99. The metallicity of the stellar popu- 
lation is taken to be solar. We considered a continuous star 



5 A short discussion of a model with more elaborate geometry is 
presented in Appendix [Bj where we show that the main conclu- 
sions arc indeed not affected. 
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Figure 3. How the photoionzation models are computed. Each 
cloud is photoionized by radiation coming from the galactic disk 
(OB stars) and from the thick disk/halo (HOLMES). The com- 
putations are done in plane parallel geometry with the ionizing 
source on one side only. 



formation and took the radiation field at an age when it is 
stabilized (in practice 10 7 yr). We do not consider dilution 
of the OB radiation field, since the latitudes we consider 
(up to 4 kpc) are much smaller that the radius of the stellar 
disk in NGC 891, which is of ~2 1 kpc (|van der Kruitlll984l ; 
iTikhonov fc Garazutdinovall2005l ). 

We estimate the total number of hydrogen ionizing pho- 
tons, Q(H°)ob, in the following way. The total infrared lu- 
mino sity of NGC 891 is Lir = 2.6 x 10 10 L (|Gao fc Solomon! 
2004). Using their Eq. 9 we obtain a total present-day star 
formation rate of 5.2 M Q yr" 1 in NGC 891 (this assumes 
that the observed FIR emission originates primarily from 
dust heated by O, B, and A starsjj. The total number of 
hydrogen ionizing photons for a constant stellar mass for- 
mation rate of 1 Mq yr _1 given by Starburst99 with the pa- 
rameters listed above is 10 53 18 ph s _1 . Thus the total num- 
ber of hydrogen ionizing photons in NGC 891 is Q(H°)ob 
= 7.8 x 10 53 ph s _1 . The vast majority of these ionizing pho- 
tons is expected to be produced by OB stars located in the 
disk. Assuming that all of them reach the eDIG, their surface 
flux at the bottom of a cloud, $ob , would then be given by: 
$ob = <2(H°)ob /(2 tyR 2 ) = 3.1 x 10 7 (21/.R) 2 ph cm" 2 s" 1 
if the OB stars are distributed uniformly in a disk whose 
radius in kiloparsec is R. The factor 2 in the denominator 
accounts for the fact that the radiation produced by the OB 
stars must feed both the "upward" and the "downward" di- 
rections. 

The H II regions located in the disk are largely ion- 
ization bounded, and photon leakage occurs through small 
"holes" , as suggested by the th in Ha filaments se en perpen- 
dicular to the disk of NGC 891 l|Rossa et al.ll2004l ). Interven- 
ing optically thick eDIG clouds further reduce the number 
of photons from OB stars reaching the cloud under consid- 



6 As mentioned by iKennicuttJ Jl99ST) . the conversion from total 
infrared luminosity to star formation rate is only approximate 
for spiral galaxies, due to complex star formation history and 
geometrical issues, but only a rough estimate is sufficient for our 
purposes. 
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oration. The range of values of $ob explored in the models 
below takes into account the, a priori ill-defined, covering 
factor. An additional factor to take into account is the par- 
tial absorption (and hardening) of the OB photons travel- 
ling in the eDIG before reaching the cloud. This is treated 
in Appendix [C] 

4.2 HOLMES 

While Starburst99 is optimized to provide a good descrip- 
tion of the ionizing radiation field from massive stars, in 
particular through its use of modern model atmospheres 
for OB stars, it does not compute the radiation field from 
post-asymptotic giant branch stars. We therefore resort 
to use the evolutionary spectral synthesis code PEGASE 
|Fioc fc Rocca-Volmerangg 1 19971 ) with which we are able 
to compute the ionizing spectral energy distribution of 
HOLMES. Taking again a Kroupa stellar initial mass func- 
tion and a solar metallicity, we consider the spectral energy 
distribution of a coeval population of stars at an age of 10 
Gyr, i.e. approximately the age of old stellar populations. 
The real star formation history of the old stellar populations 
in the halo and thick disk of NGC 891 is probably more com- 
plicated than that of the assumed instantaneous starburst. 
However the integrated radiation field from all the HOLMES 
in an old stellar population does not depend strongly on the 
details of the star formation history, as soon as ages larger 
than 10 8 yr are concerned. Roughly, the stellar energy distri- 
bution is equivalent to that of a 10° K blackbody (thus much 
harder than the radiation field from massive stars), while the 
lu minosity remains almost co nstant, as can be seen in Fig. 2 
of lCid Fernandes et al.l (|201ll ). For the scenario we adopted, 
the total number of ionizing photons of the HOLMES (ob- 
tained by integrating the spectral energy distribution from 
PEGASE) is 7 x 10 40 ph s" 1 M ^. If we consider that the 
mass of the o ld thick disk i n NGC 891 is 3 x 10 10 M Q as 
inferred from Ivan der KruitJ (| 19901 ). the total ionizing pho- 
ton number of the HOLMES present in the old thick disk is 
Q(H°)holmes = 2.1 x 10 51 ph s _1 .0 Considering the stars 
located in only half of the halo, the predicted mean surface 
flux of the ionizing photons from the HOLMES is <1>holmes 
= 8.4xl0 4 (21/i?) 2 phs^ 1 cm -2 , since the scale l ength of the 
thick disk is very close too that of the thin disk jlbata et al.l 
2009). For simplicity, we disregard the HOLMES that are 
present in the extended halo, so that, from this point of view, 
the value of $holmes that we consider is a lower limit. As 
for OB stars, we do not consider dilution of the radiation 
field. In addition, we assume that the radiation field from 
HOLMES is the same whatever the distance from the galac- 
tic plane, in particular we do not consider any metallicity 
dependence. As a matter of fact, a metallicity of one tenth 
of solar, rather than solar, could perhaps be more represen- 
tative of the metallicity of the thick disk an d halo stellar 
populations in NGC 891 (|Reikuba et alj|2009h . but the ion- 
izing radiation field from the HOLMES is virtually indepen- 
dent of metallicity, as shown in Fig. 2 of lCid Fernandes et al.l 

7 This corresponds to about 10 4 "typical" planetary nebulae nu- 
clei emitting each 2 X 10 47 ph s _1 . The real number of HOLMES 
is probably significantly larger, since many of them are expected 
to be less luminous than a "typical" planetary nebulae nucleus. 



|201ll ). We also assume that the HOLMES radiation field is 
not modified by intervening absorption. This is a simplifying 
hypothesis which is justified by the fact that HOLMES are 
present everywhere in the eDIG. We are aware that there is 
much uncertainty involved in the predictions of stellar pop- 
ulati on synthesis mode ls after the asymptotic giant phase 
(e.g. iMarigo et al]|2010l ). However, we believe that our pro- 
cedure should roughly give the correct energy distribution 
and ionizing luminosity from the integrated population of 
HOLMES. 



4.3 The grid of photoionization models 

The photoionization models for the eDIG were com- 
puted with Cloudy version 07.02.01, last described by 
iFerland et alj (|l998l ). The models are computed with the 
plane parallel approximation. Each one is defined by the 
value of <E>ob, using the unabsorbed spectral energy distri- 
bution from Starburst99, the ionization parameter U and 
the chemical composition of the gas. So far, we have no di- 
rect information on the values of the chemical composition 
of the gas in the eDIG, and there is a priori no reason to 
assume that it is solaiQ The abundances of the heavy ele- 
ments (except N, see below) relative to O are fixed to their 
solar values as implemented in Cloudy (see Cloudy's manual 
for details). The abundances of Mg, Si, Fe are depleted by 1 
dex. The N/O ratio, a priori, has no reason to be constant, 
since the nucleosynthesis sites of N and O are very different. 
Therefore, we also vary N/O in our models. This parameter 
has little impact on the [O m]/H/3 and [O n]/H/3 ratios, 
but, of course, strongly affects the computed [N n]/[0 n] 
ratio. 

The ionization parameter U is defined by $totai / (n e c), 
where $totai = $ob + ^holmes, n e is the electron density in 
the cloud, and c is the speed of light. The value of $holmes 
is the one given in Sect. 14.21 The values of U defining the 
models are varied by changing the density once $ob is fixed. 

We compute a medium resolution (MR) grid with the 
combinations of parameters given in Table [T] In total, this 
makes 34 02 models, all o f which are included in the 3MdB 
database (|Morisseul2009l ) under the reference DIG_MR. To 
improve the precision of the model-fitting for each value of 
z, we produce a high resolution (HR) grid by interpolating 
between the values obtain from the MR grid, dividing the 
size of all the steps by 2. This leads to a grid with 16 times 
more values than the MR grid, i.e. a total of 54432 models, 
referenced in 3MdB as DIGJHR. To speed-up computations, 
all the models are dust-free. The (negligible) effect of the in- 
teraction of dust with photons and gas particles is discussed 
in Appendix lAl 

Figure [4] shows the spectral energy distribution of the 
radiation emitted by the massive OB stars (red dashed 
curve) and by the HOLMES (black curves) determined as 
explained in Sects. 14.11 and 14.21 From this, one can infer 
that models with larger values of $holmes /$totai will have 
a harder ionizing radiation field and will thus produce larger 



8 Varying the metallicity of the HOLMES or OB stellar models 
would be an additional refinement, however not justified by the 
degree of approximation of our study. 
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Table 1. Parameter values for our medium resolution photoion- 
ization model grid. 



Parameter 


Min 


Max 


Step 


N a 


log $ B 


3.5 


7.5 


0.5 


9 


log O/H 


-4.3 


-2.7 


0.2 


9 


log N/O 


-1.4 


-0.2 


0.2 


7 


log U 


-4.0 


-3.0 


0.2 


6 



"Number of different values for each parameter. 
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1 2 3 4 5 
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Figure 4. The spectral energy distribution of the radiation emit- 
ted by HOLMES (black line) and massive OB stars (red dashed 
line). 

[O iii]/H/3 and [O n]/H/3 ratios. Note that the photons con- 
tributing mostly to the heating by ionization of hydrogen are 
between 1 and 3 Ryd, but, in the case of HOLMES photons 
with energies above 4 Ryd also contribute to the heating, 
due to the photoionization of He + ions, which, in the most 
ionized zones, are found in larger amounts than neutral hy- 
drogen particles, in spite of the He/H ratio being one tenth. 

5 THE DIG EXPLAINED! 

Figure [5] shows the location of the observational points in the 
[O iii]/H/3, vs [O n]/H/3 diagram, with respect to our grid 
of models. The dashed lines correspond to models with same 
values of log <E>ob (ranging from 3.5 to 6.5 by steps of 0.5), 
while the continuous lines correspond to models with same 
values of log U (equal to -4, -3.5, -3) . Each panel corresponds 
to a different value of log O/H (equal to -3.9, -3.3, -2.9). The 
value of N/O is the same in all the panels, namely -0.5 dex 
(and does not strongly affect the position of the models in 
this plot). The observational points that have the highest 
values ([O m]/H/3, [O n]/H/3) correspond to the smallest 
values of <£>ob . Since we do not know the value of the metal- 
licity in the eDIG a priori, different models can reproduce 
a given couple of observed values ([O iii]/H,9, [O Il]/H/3). 
It can be seen from these plots that the fitting models will 
not differ very much in terms of U but can have significantly 
different values of $ob, depending on O/H. There is a priori 
no way to discriminate between solutions with low and high 
metallicitjjf). 

9 Unlike in the case of classical H II regions where the degeneracy 
between high and low metallicities derived by strong line meth- 



Figure[6]is the same as Fig. [5] but in the ([O IIl]/H/3 
vs |N if] /Ha) pla ne (the famous BPT diagram from 
iBaldwin et al.|[l98lh . In all the panels, we show models with 
N/O = -0.5 dex. It is clear that the models fitting the ob- 
served data in this figure will strongly depend on the value 
of N/O. This illustrates how misleading the interpretation 
of observed spectra can be if ignoring the possible variations 
of N/O. 

For each observational point in the eDIG of NGC 891 we 
select the models of our finely meshed (HR) grid that simul- 
taneously reproduce the values of [O m]/H/3, [O n]/H/3, and 
[N ll]/Ha, within the observational uncertainties depicted in 

Fig.m 

In Fig. [7] we show the acceptable ranges of values for 
$HOLMEs/ ( I , totai, U, n c , N/O, and the mean electron tem- 
perature as a function of z. Here, n c is obtained assum- 
ing R — 21 kpc. If a smaller value of R were judged to be 
more appropriate, the values of the electron densities would 
be multiplied by (21/i?) 2 , but the other plotted quantities 
would remain the same. For example, adopting an effective 
radius of about 4.8 kpc as being mor e representative o f the 
distribution of stars in the thick disk (jlbata et al.ll2009T ) the 
values of n c would range from ~ 4 cm" 3 close to thin disk 
to ~ 0.2 cm" 3 at z = 2 kpc. This change of density is a con- 
sequence of having to keep U to the same value to obtain 
the same line ratios. We do not plot the acceptable range 
of O/H, since, for most positions, a solution can be found 
for a wide range of values of O/H. A broad range of possi- 
ble values of $HOLMEs/ < I > totai is estimated for any value of 
z but, if we restrict to models with solar metallicity (thick 
black bars in Fig. [TJl , we find that this parameter tends to 
increases with |z|, implying that the role of HOLMES in the 
ionization and heating of the eDIG clouds becomes more 
important away from the plane of the galaxy. [3 This pro- 
duces a rise in the electron temperature as seen in Fig. 03. 
The value of U shows no clear trend, while n a tends to de- 
crease with \z\. This latter result is in concordance with the 
measured behavi our of n e as function of gal actic height in 
the Milky Way l|Berkhuiisen fc Miillerl 120081 ). Notice that, 
for simplicity, we have kept the value of $holmes constant, 
while the space density of HOLMES is expected to decrease 
with \z\, so the decrease of the electron density is likely to 
be even stronger. The N/O ratio as a function of z shows an 
"M" shape, reminiscent of the observed [N n]/[0 n] ratio, 
but with larger error bars. Although a constant value of N/O 
is marginally acceptable for the entire z range sampled, we 
believe that the tendency we see for N/O could be real. As a 
matter of fact, we would expect N/O to be highest in the re- 
gion where HOLMES are important, if the diluted envelopes 
of post-asymptotic giant branch stars, which are enriched in 
nitrogen through dredge-up processes in the parent stars, 
make a significant contribution to the eDIG. 

Panels f-i of Fig. [7] compares the models that fit 
[O iil]/H/3, [O n]/H/3, and [N u]/Ra simultaneously with 

ods can be solve d using the [N ll]/[0 II] or [O m]/[N II] ratios 
llMcGaughlll994h . the nitrogen lines are of no use here, since a 
priori we know nothing about the N/O ratio. 
10 It is tempting to invert the argument and say that, since 
HOLMES are expected to play a role which is increasing with 
the distance to the galactic plane, the oxygen abundance in the 
eDIG must be rather uniform and not far from solar. 
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Figure 5. The positions of the observational points with respect to our grid of models in the [O m]/H/3 vs. [O n]/H/3 plane, for various 
values of O/H: from left to right, the value for log O/H is -3.9, -3.3, and -2.9. The dashed lines join models with same values of log <3?ob 
(equal to 3.5, 4., 4.5, 5., 5.5, 6. 6.5), while the continuous lines join models with same values of log U (equal to -4, -3.5, -3). 




Figure 6. Same as Fig.|5]but in the [O iii]/H/3 vs. [N Il]/Ha plane. 



the observations of a fe w other line ratios : He I A5876/H/3 
(obse rvational dat a form lRand et alj (12008]')) . [O l]/Hq (da ta 
from iRandl l|l998h ). [S n]/Ha (data from lOtte et"aD toll)) 
and [Ne m] A 15.6/im/[Ne 11] A12.8/xm (data from lRand et al.l 
|2008l . l2011r i). The He I A5876/H/3 ratio is reasonably well 
reproduced, given the uncertainty in the intensity of the 
weak He I A58761ine. The [S n] /Ha and [O i]/Hq ratios 
in the models does not reproduce the observed trends very 
well especially in the regions of small \z\. However, as men- 
tioned above, we cannot expect our constant-density mod- 
els to reproduce the entire range of excitations. In addi- 
tion, the [O l]line can be produced in photodissociation re- 
gions in the presence of a large proportion of very small 
grains which, as shown bv lStasiriska &i Szczer bal (|200ll ). can 
boost the gas temperature due to photoelectric heating 
and therefore allow the excitation of the [O 1] line even 
in a neutral medium. Such a scenario was not included 
in the present photoionization models. Regarding [Ne in] 
A15.6^tm/[Ne n] A12.8/im, our models overpredict this ra- 
tio. It must be recalled that to produce Ne ++ requires 
photons with energies of at least 3Ryd, i.e. much higher 
than needed to produce ++ . Photoionization model pre- 
dictions of the ionization structure of neon strongly rely on 
the number of stellar photons just above 3Ryd. For mas- 
sive OB stars, state of the art stellar atm osphere models 
differ significantly in this energy ra nge (see iMorisset et al.l 
|2004 ISimon-Diaz fc Stasinskall200ct ). For HOLMES, in ad- 
dition, the predicted spectral energy distributions in this en- 



ergy range strongly depend on the adopted abundances for 
the heavy elements, particularly of carbon. Therefore, we do 
not consider that the discrepancy between our models and 
observations as regards [Ne ill] A15.6^m/[Ne 11] A12.8/im 
disproves our models for interpreting the ionization of the 
eDIG. 



6 SUMMARY AND PROSPECTS 

There has been no agreement so far on the nature of the 
ionization sources that could explain the behaviour of the 
emission lines observed in the diffuse ionized gas (DIG) in 
late-type galaxies. Several scenarios have been proposed in 
the past, none of them being fully satisfactory. In this pa- 
per we have analyzed the galaxy NGC 891, the prototype 
and most studied of edge-on spiral galaxies, and we pro- 
posed a model that reproduces all the considered constraints 
satisfactorily. The most important challenge for the models 
considered so far was the observed increase of [O iii]/H/3, 
[O ll]/H/3, and [N ll]/Ha with increasing distance from the 
galactic plane. 

Our model considers both the population of massive 
OB stars and that of hot low-mass evolved stars (HOLMES) 
expected to be present in the galaxy. 

The ionizing photons from the massive stars located in 
the galactic plane partly escape into the thick disk and halo, 
where they are gradually absorbed by the clouds that con- 
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Figure 7. Models that fit [O mJ/H/3, [O n]/H/3, and [N Il]/Ha simultaneously. The light grey bars indicate the range of values for wich 
the models in the grid fit the observations at a given z. The dark bars show the same, but restricting to solar abundances models. The 
observed values are indicated by the big red dots. The various panels show the behaviour of various input or output parameters of the 
models: a) 'J'HOLMEs/'^'totai; b) n c (for a value of R equal to 21 kpc, see text); c) U; d) N/O ; e) the mean electron temperature; f) He I 
A5876/Ha; g) [O l]/Ha; h) [S Il]/Ha; i) [Ne III] A15.6/im/[Ne II] A12.8Atm. 



stitute the diffuse medium. The surface fraction of the gas 
clouds that are not shielded by intervening clouds becomes 
smaller with increasing distance to the galactic plane. At 
large distances from the galactic plane, HOLMES become 
dominant in terms of ionizing photons flux. Since their inte- 
grated radiation field is much harder than that of OB stars, 
the heating of the eDIG becomes more efficient. 

We have constructed a finely meshed grid of photoion- 
ization models in the framework of this scenario and, for 
each value of the distance z to the galactic plane, we have se- 
lected the ones that fit the observed values of the [O iii]/H/3, 
[O n]/H/3, and [N ll]/Ha ratios. Since we have no a priori 
knowledge of the chemical composition of the eDIG, we had 
to leave O/H and N/O as free parameters. 

We found that solutions exist for each value of z. As 
a matter of fact, the problem is not fully determined, since 



solutions can be obtained for almost any value of O/H. The 
models clearly indicate a systematic decrease of the electron 
density with increasing z. If we restrict to solar metallicity, 
we find that the models which fit the observations become 
dominated by the HOLMES as z increases. Turning the ar- 
gument around, this might be an indication that the metal- 
licity of the eDIG is roughly solar. Our models also indicate 
that N/O increases with increasing distance to the galactic 
plane, at least until |z|~ 1.5 kpc. This result is in agreement 
with what common sense would predict about the enrich- 
ment of the eDIG by HOLMES progenitors. 

In summary, we have provided a quantitative model to 
explain the behaviour of the emission lines in the eDIG of 
NGC 891 using what is known of the populations of mas- 
sive stars and HOLMES in this galaxy. This galaxy was 
chosen because it is the one with the largest set of obser- 
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vational data, allowing us to constrain the model relatively 
well. However, the model being at the same time quite natu- 
ral and robust, it is likely that it can be applied to the eDIG 
of other galaxies as well. 

Obvious next steps of this investigation would be to ob- 
tain deep 3D-spectroscopy of extraplanar gas with integral 
field units, in order to constitute a complete and reliable 
data set over the entire region of interest. This will allow us 
to map the characteristics of the stellar population using a 
stellar population synthes is code such as STARLIGHT (see 
ICid Fernandes et aij|2005h . and estimate directly the ioniz- 
ing radiation fields from the H OLMES popula ti on, in a way 
similar to what was done in IStasihska et al.l (|2008T l. This 
radiation field will then b e fed into a 3D phot oionization 
code, such as MOCASSIN (|Ercolano et alj|2003h which can 
account for the spatial distribution of the ionizing sources. 
By confronting the computed emission line ratios with those 
observed across the field of NGC 891, this should allow us 
to test in much more detail the scenario proposed here. 
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Figure Al. The values of [O in]/H/3 vs [O ll]/H/3 for grids 
of models with and without dust in the eDIG. The models 
are constructed with log O/H = —3.3, log N/H= —3.8, log 
U = —4., -3.5, and -3., and <I>ob = 4., 5., and 7. Red dotted 
lines join the dusty models while continuous black lines join the 
dust-frce models. 
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APPENDIX A: THE EFFECT OF DUST IN THE 
EDIG 

The eDIG certainly contains dust. In the models presented 
in this paper, however, dust was not explicitly included, 
mainly to speed up the computations by concentrating on 
the basic parameter of the problem, which is the spectral 
energy distribution of the ionizing radiation field. We how- 
ever depleted the abundances of metals (Mg, Si and Fe) in 
all our calculations, as men tioned in Sect. 14. 31 As shown by 
IShields fc Kennicutt! (|l995h , it is depletion which introduces 
the strongest perturbation due to dust in the optical line 
spectra. The effects linked to the physics of dust (selective 
absorption, heating and cooling by dust grains) only play a 
minor role in our problem, as illustrated with the reduced 
grid of dusty models shown in Fig. lAll Those models, joined 
by the red dashed lines, were computed taking the default 
Cloudy values of dust-to-gas ratio and dust composition, and 
cover the full range of U and $ob values. For comparison, 
we show the corresponding dust-free models joining them by 
continuous black lines. 



Figure Bl. The values of [O in]/H/3 vs [O ll]/RS for grids of 
models with different treatment of the HOLMES radiation field 
(see text). The black lines correspond to the same models as in 
Fig. lAll while the red dashed lines correspond to models in which 
the HOLMES radiation field has been separated into two compo- 
nents, ionizing a DIG cloud from "above" and from "below" (see 
text). 



APPENDIX B: THE EFFECT OF THE 
GEOMETRICAL DISTRIBUTION OF THE 
HOLMES RADIATION FIELD 

In the main body of this paper, in order to reduce the pa- 
rameter space, the transfer of radiation was computed in 
a simplified way, assuming that each cloud is ionized from 
the bottom by the radiation field from both the OB stars 
and the HOLMES. To show how the effect of a different ge- 
ometrical distribution of the HOLMES affects the results, 
we plot in Fig. [Bl] the values of [O in]/H/3 vs [O n]/H/3 
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Figure CI. The total number, Q(E), of photons above the en- 
ergy E, as a function of E, normalized to the same value Q(0), 
leaking out of photoionization models computed with log U= -3.5, 
for different values of the optical thickness at the Lyman edge, 
T13.6- The colors (black, red, green, blue, magenta, cyan, red, 
blue, magenta, cyan and green) correspond to T13 6 = 0., 1.15, 
2.30, 3.45, 4.6, 5.75, 6.90, 10., 12.9, 15.8, and 20 respectively, for 
the OB stars. The thick black line shows, for reference, the value 
of Q(E) for the HOLMES radiation field. 

for the same grid of models as shown in the central panel 
of Fig. [5] They are joined by continuous (black) lines. We 
then consider pairs of models constructed with the same val- 
ues of density and $ob, but now distributing the radiation 
of the HOLMES into 0.5$holmbs from "below", added to 
the radiation field from the OB stars, and 0.5$holmes from 
"above" . Those summed models are joined with (red) dashed 
lines. Figure IB Tl shows that all the models with the same 
value of $ob are found along the same curves, whatever 
the distribution of the HOLMES radiation. This is because 
what counts most, in the direction of increasing [O m]/H/3 
and [O 11] /H/3, is the global energy provided by the ioniz- 
ing photons. For a given observational point, however, the 
model best fitting the data will have a slightly different value 
of density than the simple models we considered in the main 
body of the paper. To obtain similar line ratios with different 
models, a similar ionization parameter U is needed, which 
is achieved by changing the gas density in the same pro- 
portions as the photon density. Therefore, the values of n 
corresponding to the observations are likely to be slightly 
different from the ones shown in Fig. [7] but their variation 
with \z\ should not change, and the main result of our pa- 
per, i.e. that the HOLMES expected to be present in the 
galaxy can explain the observed intensities of the emission 
lines above the galactic plane, remains unaffected. 



APPENDIX C: THE EFFECT OF PARTIAL 
ABSORPTION OF THE IONIZING RADIATION 
FROM THE MASSIVE OB STARS 

The dust-free case As mentioned in Sect. 14. ll the radia- 
tion field from the OB stars reaching the eDIG clouds may 
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Figure C2. The values of [O in] /H/3 vs [O 11] /H/3 for grids of 
models with different values of the partial absorption of the OB 
radiation field. The colors (black, red, blue, cyan, and green) cor- 
respond to Ti3.6 = 0., 1.15, 3.45, 5.75, and 20 respectively. The 
models have the same values of O/H, N/O, ^HOLMEs/^total an d 
U as in Fig. IA1I 

have been affected by partial absorption, especially in the 
disk. Figure IC1I shows how the values of Q(E), the total 
number of photons above the energy E, vary with E, for 
various values of T13.6. The absorption produces a harden- 
ing of the ionizing radiation. This leads to an increase of the 
electron temperature in the ionized gas and an enhancement 
of collisionaly excited lines with respect to recombination 
lines. Figure IC2I shows the effect of this absorption of the 
OB stellar radiation field on the same reduced grid of pho- 
toionization models as shown in black in Fig. lAll for some of 
the values of T13.6 considered in Fig. IC1I (and adopting the 
same colours). We see that the effect of this hardening on our 
photoionization models, while appreciable at low values of 
$HOLMEs/$totai, becomes negligible when <&HOLMEs/$totai 
exceeds 50%. This justifies our neglecting the effect of par- 
tial absorption in our high resolution grid used in Sect. [SJ 
since we are mainly interested in explaining the line ratios 
of the eDIG at high values of \z\. 

The grey dotted lines in Fig. IC2I correspond to T13.6 = 
20. For this value of T13.6, the hardness of the partially 
absorbed OB radiation field is equivalent to that of the 
HOLMES, as can be seen in Fig. lCll so that one might think 
that the predicted line ratios would be very similar to those 
of models photoionized exclusively by HOLMES. Interest- 
ingly, Fig. IC2l shows that this is not the case. This is due to 
the effect of fluorescent excitation of H/3 by stellar photons 
around the Lyman lines, which are numerous with respect 
to ionizing photons in such a circumstance, as can be seen in 
Fig. ED leading to a decrease of the [O 11] /H/3 and [O III] /H/3 
ratios due to a strong increase of H/3. The effect of fluores- 
cence is better illustrated by the top left panel of Fig. IC3I 
where the dependence of the [O ill] /H/3 and [O 11] /H/3 ratios 
with respect to T13.6 is shown. A single set of parameters 
(O/H, N/O, $ob, and U) is used. Two cases are shown: the 
red curve with square symbols corresponds to models taking 
fluorescence into account; the black curve with filled circles 
corresponds to models without fluorescence. When fluores- 
cence is not taken into account, one can reach the region 
of highest [O 11] /H/3 and [O m]/H/3 values observed in the 
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eDIG using the hardened radiation from massive OB stars 
only (for convenience, the observational points are plotted 
again in th e lower left panel of Fig. IC3j) . This was the hy- 
pothesis of ISokolowskil (| 19931 ), but we show here that the 
high values of the two line ratios are not accessible, due 
to fluorescence which increases the H/3 line intensities. No- 
tice that these models are performed using version c08.01 of 
Cloudy, which treats the fluorescence correctly. 

The dusty case Accounting for the presence of dust in 
the H uregions of the disk reduces the effect of fluorescence, 
since dust can absorb the stellar radiation below the 1 Ry- 
dberg limit. In the upper right panel of Fig. IC3I we explore 
the effect of the presence of dust in the H II regions of the 
disk, taking the default Cloudy values of dust-to-gas ratio 
and dust composition. The symbols correspond to the same 
values of Ti3.6 as in the upper left panel, but this time we 
take into account the absorption of both gas and dust. In 
all the models, fluorescence is taken into account. Each se- 
quence corresponds to different values of the ionization pa- 
rameter. The black curve marked with circles corresponds 
to log U — —3, the red one with squares to log U = —2 
and the blue one with triangles to log U = — 1. These curves 
mimic different stages of the development of an H II region, 
with the largest values of U corresponding to the youngest 
stages, where the gas is more compact and the effects of dust 
are important. As the H II region ages, it expands, the ion- 
ization parameter drops and the effects of dust become less 
important, even when keeping the same dust-to-gas mass 
ratio. In the case of a dust dominated attenuation (ie. log 
U = — 1), the eDIG models ionized by partially absorbed ra- 
diation from massive OB stars can reach the highest values 
of the observed [O ll]/H/3 and [O ill] /H/3 ratios. However in 
this case, the emission line ratios of the absorber, which are 
shown in the bottom right panel of Fig. IC3I correspond to 
a zone in the [O IIl]/H/3 vs [O Il]/H/3 diagram which is very 
far from what is observed in the disk (represented by the red 
points in the lower left panel). This means that H n regions 
with such high values of U have a negligible contribution to 
the overall emission in the disk and, therefore, that leakage 
of ionizing radiation from such H II regions into the eDIG 
must be negligible. 
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Figure C3. [O m]/H/3 vs [O ll]/H/3. Upper left: Models of the eDIG ionized by radiation from OB stars (no HOLMES) that suffered 
partial absorption by dust- free H II regions in the disk. The models have the same values of O/H, and N/O as in Fig, I All and log 
t/=-3.5. The symbols mark different values of the partial absorption (ti 3 . 6 = 0., 1.15, 2.30, 3.45, 4.6, 5.75, 6.90, 10., 12.9, 15.8, and 20 
from bottom to top). The black line with filled circles connects models where fluorescence is disabled, while the red line with square 
symbols connects models with fluorescence taken into account. Upper right: Models of the eDIG ionized by radiation from OB stars (no 
HOLMES) that suffered partial absorption by dusty H II regions in the disk. Fluorescence is taken into account. The black, red and blue 
curves correspond to absorbers with ionization parameters log U equal to -3, -2, and -1 respectively. The symbols (circles, squares and 
triangles, respectively) mark the same values of th e partial absorptio n as in the upper left panel, but this time T13.6 includes absorption 
by dust grains. Lower left: The observations from lOtte et alj i200ll ). Those corresponding to \z\ < 0.3 kpc are marked in red, the others 
in black. Lower right: Calculated emission in the disk from H II regions responsible for the absortion of the ionizing continuum used 
by the eDIG models shown in the upper right panel. The colors and symbols are the same (except the meaningless case Ti3.6=0, which 
is not shown). 



